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Soliton-pair nucleation rates F — Aexp(-B) are studied in highly-biased sine-Gordon systems 
with a local inhomogeneity for both a thermal activation regime and a quantum tunneling regime. 
It is found that the local inhomogeneity strongly affects the nucleation rates by modifying the 
bias-dependence of the exponent B. This change in the exponent B is explained as a dimensional 
crossover caused by the local inhomogeneity. It is also shown that there is another crossover, at 
which A becomes independent of the system size. 
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The dynamics of a one-dimensional sine-Gordon (SG) 
model has provided conventional understandings of vari- 
ous systems in physics, including dislocation in crystals, 
charge density waves (CDWs) in quasi-one-dimensional 
materials, and long Josephson junctions. One of typi- 
cal processes in these systems is nucleation of soliton- 
antisoliton pairs driven by an external force. The nu- 
cleation rates of soliton-pairs have been calculated in 
both a thermal activation regime and a quantum tun- 
neling regime for homogeneous systems a nun Inhomo- 
geneities, however, are unavoidable in actual experimen- 
tal situations. They may change the nucleation process 
drastically, because the nucleation rate T = Ae~ B is very 
sensitive to local modulation in the exponent B, which 
is proportional to a barrier height. Especially, inhomo- 
geneities may change the bias-dependence of the expo- 
nent B, which can be detected experimentally. This pos- 
sibility has been pointed out first in the context of CDW 
systems.© 

In this Letter, effects of a local inhomogeneity on 
soliton-pair nucleation are studied at a high driving bias. 
The nucleation rates are evaluated by Langer's methodliP 
The bias-dependence of B is typically expressed as B oc 
(fc — /) 7 j where / and f c are an external bias and a 
classical threshold bias, respectively. It is shown that 
the inhomogeneity changes not only the threshold cur- 
rent f c , but also the exponent 7. This change of 7 affects 
the nucleation rates drastically, and is expected to be de- 
tected experimentally. It is claimed that the change of 
7 is essentially understood by a dimensional crossover 
caused by the local inhomogeneity. Another crossover 
about the prefactor A is also discussed. Although the 
following discussion is available for various physical sys- 
tems described by the SG model, long Josephson junc- 
tions (LJJs) are considered as a comprehensive example, 
which allows well-controlled experiments. 

The classical equation of motion of LJJs with a local 



* Submitted to J. Phys. Soc. Jpn. 



impurity is given as 



sin <f> — f — eS(x) sin <j) — 0, 



(1) 



where / is an external current density normalized by the 
critical current density, and e(> 0) is a strength of an 
impurity potential made by modifying the thickness of 
insulator layers locally. Here, the spatial and temporal 
variables are normalized by the Josephson length and 
plasma frequency in LJJs. In this Letter, dissipation due 
to quasi-particle currents is assumed to be small, but not 
extremely weak, to guarantee the thermal equilibrium 
in a metastable state. The partition function of this 
system is described by the imaginary-time path integrals 
as Z = JVcj)(x, t) exp(— S E /g 2 ), where 



S E [<t>(x,r)] 
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is the Euclidean action and g 2 is the normalized Planck 
constant. Here, L and T are the length of the junction 
and the temperature normalized by the Josephson energy 
per unit length, respectively. At high biases / = 1 — r] 
(77 <C 1), the potential energy is allowed to be expanded 
to a quadratic-plus-cubic form. By changing the field 
variable as cf)(x) = ■K/2 + y/2r)({p{x) — l), and by rescaling 



L _„ -1/4, 

(2?7) _1 / 4 r'), the Euclidean action is modified as 
S e Mx',t')} = (2 V ) 1 I d.r' / dr' 



the spatial (temporal) variables as x = (2ry) 
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Here, e — e(2r))~ 3 / 4 is an effective impurity strength, 
which determines the magnitude of the inhomogeneity 
effects on nucleation. It should be noted that the im- 
purity effects can be controlled by the external current 
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/ = 1— 7]; the effective impurity strength e = e(2ri)~ 3 / 4 is 
enhanced if the current approaches the classical thresh- 
old current as r\ — » 0. 

The nucleation rates are evaluated by Langer's method 
in terms_oi _thii imaginary part of the free energy F — 
—T\nZQ aai3 The partition function Z is evaluated 
by integrating out the field tp(x',T') up to the second 
order of fluctuations around stationary solutions deter- 
mined by SSe/Sip — 0. In the present case, there are two 
stationary solutions: one is a stable solution ipo(x',T') 
and the other one is a bounce solution ipa{x' , r') with 
one unstable mode. The evaluated partition function Z 
includes an imaginary part produced by the integration 
around the bounce solution ip-Q^x'). The nucleation rate 
is then related to the free energy through T — 2f(T)lmF. 
Here, f(T) is a temperature-dependent factor, and takes 
1 for T < T , and T /T for T > T , where T Q is the 
crossover temperature between the thermaL activation 
regime and the quantum tunneling regime B Bill Eii The 
nucleation rate T — Aexp(-B) is then obtained as 

f(T)T — ( x(0) X 1/2 



r 



n 



(B) 



B = (S[<p B (x',T')]-S[<p (x',T')])/9 2 



(4) 
(5) 



Here, X^s (A| B ^s) are the frequencies of eigenmodes 
around ipo (<Pb) obtained by solving the 'Schrodinger' 
equation 



= X 



-dr-r- +{l~ipov{x',T'))]i>{x',T') 

(0 < B) ^v). 



(6) 



If there is a zero-frequency Goldstone mode around the 

( B ) 

bounce solution (A- — 0), this mode must be replaced 
by the_translational mode by using Fadeev-Popov tech- 
nique, a 

Thermal activation regime. At high temperatures 
T > T ~ g 2 V 1 ^ 4 i the bounce solution (pb(x',t') is 
independent of the imaginary time r', and the ImF 
method reproduces pthp-jKramers-type nucleation rate 
T = Acxp(-Af7/T)Jl3'EJ The energy barrier AU is cal- 
culated as 



AU 

(2?7) 5 / 4 



U[<pb&)] - U[y {x% 

2 



dx' 



fl 
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eipS(x') 



The stationary solutions are obtained as <Pb(x') 
ip(x';ai) and ipo(x') = ip(x';a,2) respectively, where 

3 



<p(x';a) = 



cosh 2 ((>'| +a)/2)' 



(7) 
.(8) 

(9) 



The values of a%, a2 (ai < 02) depend on the effective 
impurity strength e through e = 6a (1 — a 2 ) with a — 
tanh(a/2). At the current satisfying e = e(2ri)~ 3 / 4 — 
4/vo, the stationary solutions disappear, and the energy 
barrier AU becomes zero. Hence, the classical threshold 
current is modified by the impurity from the homoge- 
neous case rj c — as rj c — (v / 3e/4) 4/ ' 3 /2. 

The bias-dependence of the energy barrier AU ob- 
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Fig. 1. The solid line denotes the energy barrier AU as a function 
of the bias current r\ = 1 — / obtained numerically for e = 0.1. 
The dashed line denotes the result (|lo|) for e = e(2r?) -3 / 4 < 1, 
and the gray dashed line denotes that of the homogeneous case. 
The inset shows the behavior of AU near the threshold current 
V = Vc {Vc = 0.076 for e = 0.1), and the dashed line in the inset 
denotes the result fl for £ = S (e < 1). 



tained numerically for e = 0.1 is shown by solid lines 
in Fig. 1. At biases far from the threshold bias (77 3> r) c ), 
the impurity effect is weak, since the effective impurity 
strength e is small. In this case, the energy barrier AU 
is proportional to 7/ 5 / 4 as in the homogeneous case, and 
the inhomogeneity gives only a small correction of order 
of e(^C 1) as 



AU = {2rf) 



5/4 



24 



3e + 0(e 2 ) 



(10) 



This result is shown by a dashed line in Fig. 1. When 
rj approaches the threshold bias as 77 — > rj c , the effec- 
tive impurity strength e approaches the threshold value 
4/v / 3, and the bounce solution is strongly modified. As a 
result, the bias-dependence of the energy barrier is modi- 
fied. For the strong inhomogeneity e — 4/V3— £ (e <C 1), 
the energy barrier AU can be evaluated analytically as 

3/2 / \ 5/2 



AU 



f(2ifc 



5/4 { v-Vc 



-O 



(11) 



Vg v Vc j \ Vc 

This result is shown by a dashed line in the inset of 
Fig. 1. Note that the exponent of the bias-dependence 
is changed from 5/4 (eq. ©) to 3/2 (eq. @) by the 
strong inhomogeneity. This crossover about the expo- 
nent B at 77 ~ 2?7 C is expected to be detected experi- 
mentally, and is explained by a dimensional crossover as 
discussed later. 

The prefactor A is also affected by the local inhomo- 
geneity through the change in the spectrum of A^'^s. 
Generally, this change in the spectrum does not affect 
the prefactor so strongly as the change of the exponent. 
The change of the zero-frequency mode, however, may 
produce significant effects on A. In the homogeneous 
case, this zero-frequency mode around the bounce solu- 
tion (f-eix 1 ) produces the prefactor proportional to the 
system size L. In the presence of the inhomogeneity, 
this mode has a positive frequency, and modifies the In- 
dependence of A. 
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Fig. 2. The spectrum of eigenmodes around stationary solutions 
toh is shown as a function of a = tanh(a/2). For a < l/\/3, it 
denotes the spectrum around the bounce solution if-Q(x'), while 
for a > l/v3 around the stable solution ipo(x'). In addition to 
these bound states, there exists continuum spectrum at A > 1. 



The spectrum of A^°' B ^s obtained numerically from the 
Scrodinger equation (6) is shown in Fig. 2. The re- 
gion a > l/vo corresponds to the stable solution, and 
a < 1/V3 to the bounce solution. The lowest mode 

( B^ 

A-j denotes the variable along the tunneling path. The 
second mode Aj corresponds to the translational mode 
in the homogeneous case. As seen in Fig. 2, the fre- 
quency of the second mode remains small for the weak 
inhomogeneity e< 1 (a <C 1). Hence, this mode must 
be treated carefully for e <C 1 by the Fadeev-Popov tech- 
nique. Within the first-order perturbation for e, the con- 
tribution of the second mode is treated as 



A 



(B) 



(2?7) 5 / 4 24 
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dxo exp 



At/ imp /r 

cosh 2 (Cx /2) 



(12) 

where C = \/5/24. The energy modulation AC7i mp = 
3e(2ry) 5 / 4 = 3e(2rj) 1 ^ 2 by the inhomogeneity determines 
the L-dependence. The prefactor A is proportional to L 
for A[/i m p <C T (77 <C r/ cr ) as in the homogeneous case. 
However, A becomes independent of L for A£/i mp ^> T 
(77 3> T) cr ). The crossover bias is estimated as ?y cr ~ 
(T\nL/e) 2 . Note that this crossover at 77 ~ 7y cr about A 
is independent of the crossover near rj ~ 2?7 C about the 
exponent _B. 

For the strong inhomogeneity e = 4/Vo — e (e <C 1), 
the frequency of the translational mode is so large that its 
contribution to the prefactor becomes independent of L. 

(B) (Bl 

Only the frequencies of the lowest modes, \\ and A^ , 
become small compared with the characteristic frequency 
of this system (A = 1) and the other frequencies of the 
eigenmodes. Hence, only the lowest mode is relevant to 
the nucleation process. To clarify this situation, the field 
<p{x') is truncated to a one-variable problem 



tp(x') = tp(x'; a c ) + Ci(r')ipi(x'; a c ), 



(13) 



where tp(x';a) is given in 

,3 



), and ipi(x';a c ) = 
C sinh((|.x' I + o c )/2)/ cosh a (( \x' \ + a c ) /2) is the local de- 
formation mode which becomes the zero-frequency mode 
at a = a c . The normalization constant C is obtained 
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Fig. 3. The qualitative behaviors of the nucleation rate F = 
Ae~ B (a) for r? cr 2> r\ c and (b) for r\ CT r] c . The character- 
istic attempt frequency is denoted with u>q. 



from / da;' |^i | 2 = 1 as (135/(36 - 8V3)) 1/2 , and the 
critical value a c is determined by tanh(a c /2) = 1/V3- 
The coefficient C\ has dynamics in the r-direction under 
the potential 



U[ip(x')] = const. 



(14) 



Thus, the nucleation is described by the one- variable po- 
tential produced by local deformation of the field. In this 
situation, the prefactor A cannot have the L-dependence, 
because the system size L is irrelevant to the nucleation 
process caused by the local deformation of the field <f(x'). 
Note that this potential form ( |l4| ) reproduces the expres- 
sion of the energy barrier ( |TT| ) for the strong inhomogene- 
ity. 

The result is summarized in Fig. 3 for two possible 
cases. In the case of r] CI 3> rj c , the prefactor is propor- 
tional to L for the region 2rj c < rj < r] CI , and it is inde- 
pendent of L for the other current region. In the case 
of 7/ cr <C 77c, the prefactor is independent of L for any 
■q. The crossover at 77 ~ ?7 cr disappears because ?7 cr esti- 
mated for i <C 1 is not valid for the strong inhomogeneity 
77 < 2?7 C (e ~ 1), where the prefactor A is independent of 
L for any 77. 

Quantum tunneling regime. At low temperatures T < 
To ~ S^ 7 ? 1 ^ 4 ; nucleation due to quantum tunneling is 
dominant. In this regime, the stationary solutions of 
the action (3) must be calculated by solving the 1+1- 
dimensional classical field equation. However, as shown 
in the thermal activation regime, the features of the 
bias-dependence of the nucleation rate can be discussed 
by the perturbational treatment. Hence, in this Letter, 
only the limiting cases are discussed to clarify the bias- 
dependence of the nucleation rates. 

For the weak inhomogeneity e <S 1, the exponent B is 
obtained within the first-order perturbation for i as 



g 2 



~2\ 



(15) 



where sq = 31.00 and s% — 16.43. In this region, the 
exponent B is proportional to 77 1 as in the homoge- 
neous case, and the inhomogeneity effect only appears 
as a small correction. For the strong inhomogeneity 
e = 4/V3 — e (e <C 1), the system can be truncated 



Takeo Kato 



to the one- variable problem under the potential fllij) . In 
this region, the exponent B is obtained as 

15.8 



B 



,9 Z 



(16) 



Thus, the exponent of the bias-dependence is changed 
from 1 to 5/4 by the strong inhomogeneity at the 
crossover current r\ ~ 2r/ c . 

The system-size dependence of the prefactor A can 
be studied by a discussion parallel to that of the ther- 
mal activation regime. In the homogeneous case e = 0, 
there are two zero-frequency modes related to the spa- 
tial (temporal) translational symmetry of the bounce 
soliton i^g (x',t'). In the weakly- inhomogeneous case 
e <C 1, the frequency of the temporal translational mode 
remains zero, while the frequency of the spatial trans- 
lational mode is lifted. The Fadeev-Popov technique is 
applied to this spatial mode as 
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A 



const. / dxoexp 



Mi 



(17) 



Here, the function f(xo) = J dr'ip^ (x ,t') behaves as 
f(xo) — > for |xo| — ► oo, and has a maximum value 4.784 
at xo = 0. It is found from (|l7| ) that the factor er\j g 2 de- 
termines the system-size dependence of the prefactor A. 
For er\jg 2 <C 1 (f] <C T) CT ), the prefactor is proportional 
to L, while for iij/g 2 » 1 (77 3> ?7cr)> the prefactor be- 
comes independent of L. The crossover bias is estimated 
as ?7 cr ~ (g 2 lnL/e) 4 . 

Discussion. The bias-dependence of the exponent 2? 
is summarized in Table I, where the results for small 
Josephson junctions (SJJs) are also shown for compar- 
ison. At high temperatures, the exponent B behaves 
for SJJs as B oc 17 3 / 2 , and for homogeneous LJJs as 
B oc 77 s / 4 . This difference between SJJs and LJJs comes 
from the character of the bounce solution tpB '■ the bounce 
solution has spatial dependence in the x'-direction for 
LJJs, while it is independent of x' for SJJs. In other 
words, there is a spatial dimensional crossover between 
SJJs and LJJs. Then, the bias-dependence of B is ex- 
pressed as B tx rj 3 / 2 ~ d / 4 by a spatial dimension d which 
takes d = 1 for LJJs and d = for SJJs. There also 
exists a 'temporal' dimensional crossover between the 
thermal activation regime and the quantum tunneling 
regime: the bounce solution ips has a temporal depen- 
dence in the r '-direction at low temperatures, while it is 
independent of t' at high temperatures. As a result, the 
bias-dependence of B at low temperatures is expressed 
as B oc j ? 3 /2-(d+i)/4 ) w here (d + 1) denotes the total 
dimension of the system including the r-direction. 

The inhomogeneity effects on the bias-dependence of B 
can be understood in terms of a dimensional crossover. 
The inhomogeneity modifies both the bounce solution 
ipsix' ,t') and the stable solution ipo(x',T'). Although 
these two solutions have a spatial dependence, the differ- 
ence between these solutions becomes small in the pres- 
ence of the strong inhomogeneity. As a result, the nuclc- 
ation process is effectively described by local deformation 
of the field, and the spatial dimension d is effectively re- 
duced from d = 1 to d = by the strong inhomogeneity. 
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Table I. The bias-dependence of the exponent B in three sys- 
tems: homogeneous (or weakly-inhomogeneous) long Josephson 
junction (LJJ), small Josephson junction (SJJ), and LJJ with 
the strong inhomogeneity e ~ 1. 



Therefore, in strongly-inhomogeneous LJJs, the expo- 
nent of rj in B is the same as that of SJJs. (See Table I.) 

In this Letter, soliton-pair nucleation rates F = 
Aexp(-B) have been studied for highly-biased sine- 
Gordon (SG) systems with a local inhomogeneity. It 
is found that the effective inhomogeneity strength e ~ 
e(l — /)~ 3 / 4 controls the bias-dependence of the exponent 
B which is typically written as B cx (/ c — f) 3 / 2 ~ d / 4 in the 
thermal activation regime and as B cx (/ c — /) 3 / 2 ~( d+1 )/ 4 
in the quantum tunneling regime, where f c and / are 
the classical threshold bias and the external bias, re- 
spectively. The spatial dimension d takes d = 1 for the 
homogeneous or weakly-inhomogeneous one-dimensional 
SG systems (e 1), while it is reduced to d = 
for strongly-inhomogeneous SG systems (e ~ 1). This 
change in B is. expected to be detected experimentally, 
e.g. in LJJs.lisEJj) This phenomena would be available 
to evaluate inhomogeneity strength in real experimental 
systems. It is also found that there exists a different 
crossover bias rj CT at which the prefactor becomes inde- 
pendent of the system size L. These results are universal 
for systems where the nucleation process is most domi- 
nant. Details of the relevance to actual experiments will 
be presented elsewhere. 

The author thanks M. Yumoto for stimulating discus- 
sion. The author is supported by Research Fellowship 
of Japan Society for the Promotion of Science for Young 
Scientists. This study was supported by a Grant-in-Aid 
for Scientific Research from the Japanese Ministry of Ed- 
ucation, Science, Sports and Culture. 



[1] K. Hida and U. Eckern: Phys. Rev. B 30 (1984) 4096. 
[2] S. Nakaya and K. Hida: J. Phys. Soc. Jpn. 55 (1986) 3768. 
[3] B. I. Ivlev and V. I. Mel'nikov: Phys. Rev. B 36 (1987) 6889. 
[4] P. Hanggi, F. Marchesoni and P. Riseborough: Europhys. 

Lett. 13 (1990) 217. 
[5] M. Yumoto, H. Fukuyama, H. Matsukawa and N. Nagaosa: 

the proceedings of LT22, to be published in Physica B (2000). 
[6] D. W. McLaughlin and A. C. Scott: Phys. Rev. A 18 (1978) 

1652. 

[7] J. S. Langer: Ann. Phys. (N. Y.) 41 (1967) 108. 
[8] C. G. Callan and S. Coleman: Phys. Rev. D 16 (1977) 1762. 
[9] I. Affleck: Phys. Rev. Lett. 46 (1981) 388. 
[10] U. Weiss: Quantum Dissipative Systems (World Scientific, 

Singapore, 1993). 
[11] H. Grabert and U. Weiss: Phys. Rev. Lett. 53 (1984) 1787. 
[12] A. I. Larkin and Yu N. Ovchinnikov: Zh. Eksp. Theor. Fiz 86 

(1984) 719; Sov. Phys.-JETP 59 (1984) 420. 
[13] H. A. Kramers: Physica (Utrecht) 7 (1940) 284. 



Dimensional Crossover by a Local Inhomogencity in Soliton-Pair Nucleation 



5 



[14] P. Hanggi, P. Talkner and M. Borkovec: Rev. Mod. Phys. 62 
(1990) 251. 

[15] A. Davidson, B. Dueholm, B. Kryger and N. F. Pedersen: 

Phys. Rev. Lett. 55 (1985) 2059. 
[16] A. V. Ustinov, T. Dodercr, R. P. Huebener, N. F. Pedersen, 

B. Mayer and V. A. Oboznov: Phys. Rev. Lett. 69 (1992) 

1815. 



